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SUMMARY 

This  report  develops  the  theory  for  a  multitarget  tracking 
algorithm  based  on  Probabilistic  Data  Association  with  new 
selection  rules  for  assigning  sensor  measurements  to 
target  tracks  and  for  forming  multitrack  clusters.  These 
new  rules  remove  the  requirement  to  form  a  gate  about  each 
target's  predicted  position  for  the  selection  of  sensor 
measurements .  The  resultant  algorithm  is  the  same  for  all 
target  tracks  and  clutter  conditions.  The  algorithm 
adapts  to  the  sensor  measurements  via  probability  terms 
which  model  the  environment  and  sensor  processing. 
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1 .  INTRODUCTION 

A  solution  to  the  problem  of  tracking  a  known  object  in  a  cluttered 
environment  was  addressed  by  Bar-Shalom  and  Tse(ref.2)  who  formulated  a 
sub-optimal  solution  referred  to  as  the  probabilistic  data  association  (PDA) 
filter.  Their  solution  builds  on  the  Kalman  filter  by  taking  into  account  the 
event  that  sensor  measurements  used  for  updating  the  filter  could  have 
originated  from  clutter  and  noise.  In  their  analysis,  and  that  of  extensions 
to  the  PDA  filter  in  references  3,  4  and  5,  the  sensor  measurements  for 
updating  the  filter  are  all  those  located  inside  a  'gate'  which  is  centered  on 
the  target's  predicted  position.  The  differences  between  the  predicted  target 
position  and  these  measurements  are  weighted  by  the  probability  of  their 
having  originated  from  the  target.  This  leads  to  a  weighted  difference  which 
is  used  to  form  the  filtered  estimate.  This  approach  eliminates  the  need  for 
complex  testing  procedures  before  the  tracking  step  to  resolve  the  target 
origins  of  measurements . 

With  the  use  of  a  gate,  the  number  of  measurements  used  to  update  a  PDA  filter 
can  vary  dramatically  with  time  and  clutter  conditions.  In  the  case  of  the 
Joint  PDA  filter (ref .3,4)  for  tracking  multiple  targets,  the  measurements 
within  all  the  gates  of  neighbouring  tracks  are  used  for  the  probability 
computations.  Neighbouring  tracks  are  defined  as  those  whose  gates  overlap 
and  this  group  of  tracks  is  referred  to  as  a  track  cluster.  In  this  case  the 
probabilities  of  all  origins  of  measurements  from  targets  in  the  track  cluster 
are  computed.  This  results  in  a  rapid  growth  of  the  number  of  possible 
origins  with  cluster  size  and  number  of  measurements.  Because  of  computer 
limitations,  any  implementation  of  these  algorithms  has  to  include  upper 
bounds  on  the  number  of  measurements  and  the  cluster  size. 

To  remove  the  above  variability  and  thereby  derive  a  filter  algorithm  which  is 
independent  of  target  and  clutter  conditions,  new  rules,  which  were 
foreshadowed  in  references  6  and  7,  are  introduced  for  selecting  sensor 
measurements  and  for  forming  target  clusters.  In  references  6  and  7  the  gate 
selection  approach  for  the  PDA  filter  was  replaced  by  an  approach  which 
involved  selecting  the  'm'  nearest  measurements  with  'm'  fixed.  In  this 
paper,  the  analysis  in  reference  7  is  extended  from  the  single  target  case  to 
include  track  clusters  of  fixed  size.  In  this  case  a  cluster  consists  of  a 
group  of  1 t'  tracks  which  may  not  necessarily  be  'close'  to  one  another.  In 
addition,  this  paper  also  develops  a  sensor  model  and  refines  the  analysis 
contained  in  reference  7.  The  incorporation  of  the  event  to  account  for  track 
initiation  and  termination  (see  reference  7)  is  not  included  in  this  paper  but 
will  be  treated  separately.  The  emphasis  of  this  paper  is  on  the  theoretical 
development  of  a  PDA  filter  with  these  new  selection  rules. 

The  theoretical  development  is  organized  so  that  the  first  section  gives  an 
outline  of  the  target,  sensor  and  update  models.  This  is  followed  by  the 
definition  of  the  new  selection  and  clustering  rules.  The  feasible  events  are 
then  defined  and  this  is  followed  by  the  formulation  of  the  new  algorithm 
based  on  these  rules.  The  assumptions  used  in  the  development  of  the 
algorithm  are  numbered  for  ease  of  reference  and  to  help  emphasise  the 
assumptions  on  which  the  analysis  is  based.  As  reference  is  made  to  Kalman 
filter  theory  which  uses  the  symbol  'P'  for  covariance,  the  notation  P(.)  is 
used  to  denote  a  covariance  matrix,  P  and  Pr{ . }  denote  probability.  For 
continuous  variables  p(.)  denotes  probability  density  while  for  discrete 
variables,  p(.)  represents  the  probability  mass  function.  In  order  to  avoid 
unnecessary  complication  of  the  equations,  the  time  index  is  only  given  when  a 
variable  or  event  is  first  defined.  Thereafter  the  time  index  is  only 
retained  when  required  for  clarity.  In  the  case  of  scalars  which  are  derived 
from  time  dependent  vectors  and  other  scalars,  the  time  index  is  not  included 
as  this  is  implied. 
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2.  OUTLINE  OF  TARGET,  SENSOR  AND  UPDATE  MODEL 

The  model  used  to  represent  the  environment,  target  and  sensor  is  illustrated 
in  figure  1.  Here  the  sensor  receiver  derives  its  input  from  clutter,  noise 
and  targets  distributed  throughout  the  sensor's  surveillance  region.  This 
input  can,  in  the  case  of  a  radar,  be  electromagnetic  waves  generated  from  the 
illumination  of  the  sensor  surveillance  region  with  an  appropriately  modulated 
transmitter . 

When  propagation  time  between  objects  in  the  surveillance  region  and  the 
sensor  receiver  is  neglected,  the  receiver  is  equivalent  to  a  bank  of  parallel 
receivers.  Each  of  these  receivers  gives  the  output  from  one  of  the  sensor's 
resolution  cells  in  the  surveillance  region.  The  receiver  number  or  index  ' i ' 
corresponds  to  the  resolution  cell  number.  The  output  from  all  of  the 
receivers  at  time  'k' ,  then  consists  of  N  values  which  are  represented  by  the 
set  {^(k)  , .  .  ,£N(k)} ,  where  N  is  the  number  of  sensor  resolution  cells. 

The  detection  process  following  the  receiver,  searches  the  receiver  output  to 
locate  likely  target  signals.  This  process  could,  for  example,  involve 
comparing  the  receiver  outputs  against  a  threshold  and  those  C's  which  pass 
the  threshold  are  then  used  for  parameter  estimation.  This  latter  process 
normally  involves  selecting  each  Z.  from  the  previous  test  and  then  comparing 
it  with  S's  in  adjacent  resolution  cells  to  estimate  the  surveillance  region 
coordinates  of  the  source  of  the  received  signal.  Each  estimate  from  the 
detection  process  is  referred  to  as  a  sensor  measurement  which  is  represented 
by  the  vector  z^Ck) .  This  vector  contains  components  such  as  range,  bearing, 

etc  and  ' i '  is  an  arbitrary  measurement  subscript.  The  measurements  output  by 
the  detection  process  are  represented  by  the  set  {z^(k) , . . ,z^i (k)}  where  N'  is 

the  number  of  measurement  vectors  from  the  detection  process  at  time  'k' . 

The  clutter  map  process  which  follows  the  detection  process  selects  those 
measurements  in  the  clutter  region  which  is  a  subset  of  the  surveillance 

region.  The  region  can  be  based  on  the  predicted  position  of  a  track  or  on 

the  distribution  of  the  clutter  and  noise  throughout  the  surveillance  region. 
For  the  latter  case,  the  surveillance  region  is  partitioned  into  regions  of 
approximately  constant  clutter  and  noise  density.  For  either  case,  the 
clutter  region  is  made  so  that  the  probability  density  function  for  clutter 
and  noise  is  invariant.  Measurements  derived  from  clutter  and  noise  are 
collectively  referred  to  as  clutter  measurements  and  they  satisfy  the 
following  assumption. 

ASSUMPTION  NO. 1 :  Clutter  measurements  are  assumed  to  be  independent  and 
uniformly  distributed  within  a  clutter  region  which  may  or  may  not 

contain  targets .  The  total  number  of  measurements  in  V is  known  and 

represented  by  n(k)  ,  where  n(k)»l. 

The  following  analysis  is  based  on  targets  in  an  arbitrary  clutter  region  Vc 

with  volume  v  .  Extension  to  other  clutter  regions  is  by  identical  analysis, 
c 

The  X  target's  state  at  time  index  'k'  is  given  by  x^(k).  The  model  used 

for  the  sensor  measurement  step  and  the  equations  for  the  target's  state  as  a 
function  of  the  time  index  'k'  are  represented  by: 
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ASSUMPTION  NO. 2 (a):  A  sensor  measurement  z(k)  from  target  X  is  related  to 
the  target's  state  (k)  by 


z(k)  =  H(k)x^(k)  +  v(k) 


(1) 


where  H(k)  is  the  measurement  matrix  and  v(k)  is  the  measurement  noise 
introduced  by  the  above  parameter  estimation  step  in  conjunction  with  the 
sensor  parameters.  This  noise  is  represented  by  a  zero  mean  and 
covariance  R(k) . 

ASSUMPTION  NO .2(b):  The  target  state  evolves  with  time  by 


x^(k+l)  =  F(k)xx(k)  +  u(k) 


(2) 


where  F(k)  is  the  state  transition  matrix  and  u(k)  is  the  process  noise 
with  zero  mean  and  covariance  Q(k).  The  process  noise  accounts  for  target 
manoeuvres  and  model  errors. 

ASSUMPTION  NO .2(c):  The  measurement  matrix  H,  measurement  error 
covariance  R,  target  state  transition  matrix  F  and  process  noise 
covariance  Q  are  known. 

Note  that  from  equation  (1),  the  measurement  vector  does  not  contain  a  target 
subscript  since  the  measurement  origin  is  unknown.  In  this  paper,  multiple 
targets  are  considered  and  it  is  assumed  that  they  satisfy  the  following 
assumption. 

ASSUMPTION  NO. 3:  A  target's  state  is  independent  of  all  other  targets  and 
sensor  measurements  which  originate  from  targets  are  independent  of  each 
other  and  of  n. 

Independent  target  states  are  based  on  minimum  a  priori  knowledge.  As  a 
general  rule  this  condition  is  satisfied,  eg  crossing  aircraft.  However  in 
the  case  of  aircraft  flying  in  formation,  their  states  are  no  longer 
independent.  In  this  analysis  no  use  is  made  of  any  correlation  between 
target  states.  Independent  sensor  measurements  normally  apply  unless  the 
detection  process  is  overloaded  or  the  receiver's  output  from  any  resolution 
cell  contains  a  significant  contribution  from  more  than  one  target.  This 
latter  case  occurs  when  targets  are  of  the  order  of  a  resolution  cell  apart 
and  leads  to  a  bias  in  the  parameter  estimation  step.  The  closely  spaced 
target  case  can  also  lead  to  a  merged  measurement.  This  has  been  analysed  for 
the  PDA  filter  in  reference  4.  In  this  analysis,  the  merged  measurement  is 
not  considered. 

til 

The  estimate  of  the  X^  target's  state  at  time  k,  given  data  up  to  time  j,  is 
denoted  by  x.  (k|j).  The  covariance  associated  with  this  estimate  is  denoted 

A 

by  Px(k|j).  At  time  'k'  it  is  assumed  that, 

ASSUMPTION  N0.4:  For  every  target  in  the  region  V  ,  there  is  a 
corresponding  track  which  contains  a  state  estimate  x^Ckjj)  and  covariance 
Px(k|j). 
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For  the  special  case  when  there  is  no  uncertainty  in  the  origin  of  sensor 
measurements  used  for  updating  a  target's  state  estimate,  the  associated 
covariance  is  given  by  P^CkJk).  If  Assumption  No. 2  is  satisfied,  the 

discrete-time  Kalman  filter(ref . 1)  applies  and  the  updated  state  is  given  by, 


*i\(klk>  =  x^(k|k-l)  +  W^(k)yiX(k)  (3) 

P\(k|k)  =  [I  -  Wx(k)H]Px(k|k-l)  (4) 

Where  the  innovation  vector,  innovation  covariance,  gain  matrix  and  predicted 
covariance  are  respectively, 


/N/ 

yiX(k) 

=  z.(k)  -  Hxx(k|k-1) 

(5) 

sx00 

=  HPx(k|k-l)HT  +  R 

(6) 

Wx(k) 

=  Px(k|k-l)HTSx_1(k) 

(7) 

Px(k|k-1) 

=  FPx(k-l|k-l)FT  +  Q 

(8) 

Here  H  and  F  denote  H  transpose  and  F  transpose  respectively. 

When  applying  a  Kalman  filter  to  sensor  measurements  with  unknown  origins  from 
a  multitarget  and  cluttered  environment,  the  measurements  used  for  updating 

the  X  filter  are  often  chosen  from  those  contained  within  a  selection  gate. 
This  gate  is  scaled  from  the  innovation  covariance  S. (k) ,  and  placed  about  the 
/\  ^ 

predicted  state  Hxx(k|k-1).  If  more  than  one  measurement  is  contained  within 

the  gate  or  if  the  gates  from  two  or  more  tracks  overlap,  the  measurements  are 
often  assigned  to  tracks  on  the  basis  of  nearest  neighbour.  Figure  2 
illustrates  the  case  of  a  two-dimensional  region  containing  two  targets  with 
overlapping  selection  gates  and  measurements  from  clutter  and  targets.  If 
measurements  are  assigned  using  a  nearest  neighbour  rule,  then  one  possible 
assignment  is  measurement  3  to  track  1  and  measurement  5  to  track  2.  Other 
assignments  are  possible  but,  regardless  of  the  assignment,  the  filter  update 
(see  equations  (3)  and  (4))  assumes  that  the  assignment  is  correct.  To  reduce 
the  occurrence  of  non-target  measurements  updating  the  filter  and  thereby 
minimise  errors  in  the  target  state  estimate,  the  sensor  clutter  measurement 
rate  has  to  be  maintained  at  a  low  level.  In  many  applications,  maximum 
sensitivity  is  required  and  therefore  the  clutter  measurement  rate  is  often 
adjusted  to  the  maximum  rate  that  can  be  tolerated  by  the  tracking  system. 
The  following  sections  deal  with  an  alternative  approach  to  multitarget 
tracking  in  cluttered  conditions  based  on  a  PDA  filter  which  takes  account  of 
the  clutter  and  other  target  origins  for  sensor  measurements  which  update  the 
filter. 
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3.  NEW  SELECTION  AND  CLUSTERING  RULES 

The  new  sensor  measurement  selection  rules  and  track  clustering  rules 
introduced  in  this  paper  replace  the  'gate'  concept.  These  rules  are  based  on 
a  nearest  neighbours  approach  where  more  than  one  neighbour  is  selected.  Here 
the  'm'  nearest  measurements  and  the  ' t  —  1 ’  nearest  tracks  to  each  reference 
track  are  selected  with  'm'  and  't'  constant.  In  the  following  development,  a 
track  is  considered  with  respect  to  the  ' t-l'  surrounding  tracks  and  the  track 
subscript  is  re-ordered  so  that  \=1  represents  the  reference  track.  To 
simplify  notation,  the  track  subscript  is  omitted  when  dealing  with  the 
track.  The  rules  for  selecting  the  measurements  and  the  track 
clusters  for  each  reference  track  can  be  written  as; 

(a)  Select  the  next  track  as  reference  and  set  its  track  subscript  X=l, 

(b)  Select  the  't-l'  nearest  target  tracks  based  on  the  distance, 

~T  -1~  (q\ 

TX  =  X  Xs  Xx  C9) 

where  xx  =  Hxx(k|k-1)  -  Hx^kjk-1).  The  track  subscript  ’  X ’  is  re-ordered 
so  that  X=t  is  the  most  distant  track. 

(c)  Select  the  'm'  nearest  measurements  to  the  reference  track  based  on 
the  distance. 


d. 

l 


The  measurement  subscript  ' i ' 
the  track  subscript  T  X  ’  in  equation  (10),  which  for  this  case  is  1  ,  has 
been  omitted. 

Thus  'm*  measurements  are  always  selected  for  each  track  and  all  tracks  are ^ in 
a  cluster  of  't'  tracks.  The  following  analysis  is  based  on  t<m<n.  The  'm' 
nearest  measurements  to  the  reference  track  represent  an  ordered  set  derived 
from  the  the  set  of  'n'  measurements  in  the  region  V^.  The  ordered 

measurements  are  denoted  by  y^(k) ,  i=l,..,m  with  the  ordered  set  defined  by, 

Y(k)  =  {y1(k),..,ym(k);  d^.  .<dm>  (11) 

For  each  track  in  the  cluster,  the  past  data  up  to  time  k-1  are  denoted  by 

T  (k-1).  All  the  track  data  in  the  cluster  of  't'  tracks  are  represented  by 
X 

the  set. 


=  yT.S_1y.  (10) 

J  i  3  i 

is  re-ordered  so  that  i=l  is  the  nearest  and 


•  ,t 


T(k-l) 


(12) 
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From  these  definitions,  the  same  set  of  'm'  measurements  are  not  necessarily 
selected  when  the  other  't-l'  tracks  in  the  cluster  are  used  as  the  reference 
track.  This  is  not  considered  to  result  in  any  penalty  because  the  selected 
measurements  form  the  group  of  most  probable  measurements  for  each  reference 
track.  Also  this  new  procedure  requires  there  to  be  ' t '  tracks,  ie, 

ASSUMPTION  NO. 5:  There  are  at  least  '  t'  targets  in  the  clutter  region  V  . 

Measurements  from  targets  which  are  not  included  in  any  cluster  of  1 t' 
tracks  are  equivalent  to  clutter  measurements. 

To  ensure  that  there  are  always  ' t '  tracks,  dummy  tracks  can  be  included  with 
parameters  such  that  they  do  not  interfere. 


4.  EVENT  SPACE  FOR  SELECTED  MEASUREMENTS 

The  measured  values  at  time  'k'  from  each  of  the  process  outputs  in  the 
previous  section  represent  the  values  of  random  variables.  Throughout  this 
section,  a  symbol  without  a  time  index,  eg  £^,  denotes  a  random  variable  while 

a  symbol  with  a  time  index  is  the  value  of  the  random  variable.  Thus  £^(k)  is 

the  value  of  the  random  variable  £^  at  time  'k' .  Now  the  N  receiver  outputs 

from  the  surveillance  region  represent  a  random  sample  {£^,..,£^}.  Each  £ 

originates  from  either  clutter  or  target  plus  clutter  in  the  sensor's 

resolution  cell.  To  denote  the  origin  of  the  i^  receiver  output,  an  origin 
index  2^  is  introduced.  This  origin  index  has  the  range  0,1,.., t  with  2^=0 

used  to  denote  a  clutter  origin  and  with  £>0  used  to  denote  a  target  plus 

clutter  origin  associated  with  track  X=2^.  A  target  plus  clutter  origin  is 

th 

referred  to  as  a  target  origin.  Therefore  the  pair  denotes  the  i 

receiver  output  with  origin  defined  by  £^.  Because  the  numbers  of  the 

non-zero  origins  are  identical  to  the  re-ordered  track  subscripts,  an  origin 
can  be  used  for  a  track  subscript. 

The  elements  r  of  the  receiver  output  space  R  consist  of  N  receiver  outputs 
with  N  origins  and  past  target  data  T.  By  grouping  the  £  and  2  variables 
separately,  the  element  r  is  given  by  r=[ (£^, . . ,£^) , (2^, . . , 2^) ,T] .  The  output 

space  R  is  then  constructed  from  all  possible  values  for  the  element  r.  These 
values  are  obtained  by  allowing  T  and  the  N  variables  £  and  2  to  vary  through 
all  permissable  values.  Remember  that  the  origin  value  of  each  component  in 
the  elements  of  the  sample  space  is  known.  The  valid  lists  of  origins 
contained  in  each  element  r  have  a  maximum  of  one  occurrence  of  each  target 
plus  clutter  origin.  Now  the  measurements  with  a  clutter  origin  are  samples 
from  a  common  distribution  function  whereas  each  measurement  with  a  target 
origin  is  from  a  unique  distribution  function  which  is  indexed  by  the  origin 
value.  From  the  definition  of  the  event  space,  there  are  elements  with  the 
values  £  (k),..,£  (k)  and  associated  origins  2^ (k) , . . , 2^(k)  which  have 

identical  £  and  2  values  but  the  locations  of  the  clutter  measurements 

re-arranged.  These  elements  have  the  same  joint  distribution  function  on  R. 

No  similar  re-arrangement  of  target  origin  subscripts  exists  because  each 

element  r  has  a  maximum  of  one  measurement  for  each  2.>0.  Remember  that  each 

i 

arrangement  of  clutter  subscripts  represents  a  different  element  r. 
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The  detection  process  which  follows  the  receiver  (see  figure  1)  is  represented 
as  a  mapping  from  the  receiver  output  space  R  to  the  detection  output  space  U. 
A  detection  space  element  is  given  by  u=[  (z.^ ,  .  .  .z^,  )  ,  (^ ,  .  .  ,  )  ,T]  where  the 

range  of  N'  is  m+l,..,N.  The  mapping  is  illustrated  in  figure  3  where  the 
vector  function  h(r)  is  used  to  denote  the  detection  process.  From  Section  2, 
the  detection  process  operates  on  a  number  of  receiver  outputs  for  the 
parameter  estimation  step.  Therefore  each  measurement  with  origin  attached  is 
given  by  (z . , £ . )=h . ({?, 1}  .  )  where  £.=£.  and  {?,£}.  is  the  set  of  receiver 

outputs  used  for  parameter  estimation.  The  components  of  the  elements  of  the 
sample  space  U  are  obtained  from  those  receiver  outputs  where  hU?,*}^. 

When  a  receiver  output  is  obtained  with  a  target  origin,  the  target  is 
detected.  Detection  is  represented  by  the  event  D£  which  is  defined  by, 


=  (u=h(r) ;  (z . ,£i=£)=hi(r)  for  some  i=l,.,N',  and  VrtR}  £-0,1,.., t  (13) 


Where  h.(r)  is  the  it^1  component  of  the  vector  function  h(r) .  From 
equation1  (13),  D£  and  D £ ,  for  each  £>0,  partition  U  and  figure  3  gives  the 
case  for  one  target  only.  In  the  case  of  DQ,  every  element  of  U  has  at  least 
one  component  with  a  clutter  origin.  Therefore  DQ=U.  In  this  analysis  there 
are  '  t  *  targets  and  their  events  D£  are  independent  of  each  other  (see 
Assumption  No. 3).  For  an  arbitrary  element  ueU,  then  ueD£  or  D£  for  each  £>0. 
All  the  valid  combinations  of  D£  or  D£  for  the  V  targets  can  be  derived  from 
a  binary  number  of  V  bits.  This  binary  number  is  represented  by  the  set 


B  (t) 

K 


=  ( (b. 


’V 


bie(0,l) 


and  ic=l+Ib.2Cl'1)} 
1 


k=1,..,21 


(14) 


By  equating  the  bit  b  subscript  '  i*  to  the  event  D  subscript  l  , 
values  then  denote  D  and  D  where  b_.=l  denotes  D..  and  b.,-0  denotes  D... 

case  of  t=2,  there  are  4  partitions  on  U  where  each  partition 
intersection  of 


the  bit 
For  the 

is  the 


K 

V2) 

Event  on  U 

1 

0  0 

®1  11  ®2 

2 

1  0 

D1  "  B2 

3 

0  1 

Df  n  d2 

4 

1  1 

D1  n  D2 

Note  that  through  the  action  of  the  clutter  map  process  and  the  selection 
process,  these  partitions  map  onto  the  output  sample  spaces  of  these 
processes.  The  elements  contained  in  the  mapped  partitions  are  not  the  same 
but,  to  avoid  a  burgeoning  notation  problem,  the  same  symbol  is  used.  The 
sample  space  associated  with  a  mapped  partition  is  apparent  from  its  context. 
Furthermore,  the  origin  accompanying  the  variable  from  each  process  is 
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represented  by  the  same  symbol  £.  It  is  to  be  remembered  that  the  origin 
value  SL(k)  linked  with  say,  the  receiver  output  value  ^(k),  is  not  the  same 

origin  for  the  component  i  from  say,  the  detection  process. 

The  clutter  map  process  which  follows  the  detection  process,  selects  those 
measurements  in  the  region  V  and  is  represented  by  the  function  u(z  .,£.). 

C  11 

This  function  maps  from  0  to  Z  as  illustrated  in  figure  3.  The  elements  of 

the  clutter  map  output  space  Z  are  given  by  z={(z'  ,..,z'  ),(£,,..,£  ) ,T} 

lnln 

where  n  has  the  range  m+l,..,N'  and  z=u(u).  As  in  the  detection  process,  the 
clutter  map  does  not  pass  all  measurements  and  the  measurement  subscripts  in 
any  element  z  are  arbitrarily  assigned.  When  a  measurement  with  origin  £  is 
output,  the  event  occurs  and  is  defined  as 

V  =  (z=d(u);  (z' . ,£ ,=£)=u . (u)  for  some  i=l,.,n, 

~  111 

andVueD^}  £=0,1,.., t  (15) 

Where  u  .  (u)  is  the  itd  component  of  the  vector  function  \)(u)  .  From  the  above, 
for  origin  £> 0,  each  element  zeD^  fl  has  £=£^  for  some  i.  The  compliment 
also  exists  which  contains  those  elements  z  where  the  target  origin  £#£^  for 
all  i,  ie  zsD^flV^.  Therefore  the  binary  number  B^(t)  can  also  be  used  to 
obtain  all  possible  combinations  of  these  events  on  Z. 

The  clutter  map  output  space  Z  is  the  input  to  the  selection  process  which 
selects  the  'm'  nearest  z'  values  using  the  scalar  distance  d  in 
equation  (10).  The  selection  process  which  is  represented  by  the  vector 
mapping  function  ^(z) ,  involves  a  partial  ordering  of  Z  into  elements  with  the 
'm'  smallest  d  values.  The  remaining  n-m  z'  values  have  d  values  greater  than 
the  'm'  smallest  and  are  discarded  while  the  value  of  n  is  retained.  Thus  the 
selection  process  output  space  Y  has  the  elements 

y=[((y15..,ym), V’n>T);  di£-^dJ- 

For  any  element  zzZ,  origin  lists  'L',  'A'  and  'A'  are  introduced.  The  origin 
list  L  contains  the  origins  associated  with  the  measurements  which  have  the 
'm'  smallest  d  values.  The  origin  list  A  contains  target  origins  not  in  L  but 
contained  in  the  element  z.  The  origin  list  A  contains  the  target  origins  in 
LU  A.  In  addition  to  these  lists,  an  arrangement  list  'A',  which  contains  the 
measurement  subscripts  associated  with  the  origins  in  the  list  L  is  included. 
Thus  A  provides  the  link  between  the  origins  in  L  and  the  measurements  in  z. 
The  lists  L,  A,  A  and  A  are  then  defined  as 


L  =  {(£  ,..,£  );  d  <..<d  } 

1  m  1  m 


A 


{(x 


1*  ' 


•  ,x  );  d 

m  i 


<..<d  } 

1  lm 


A  =  the  unordered  list  of  target  origins  in  Lfl{£  ,..,£  } 


A  = 


the  unordered  list  of  target  origins  in 


LU  A 


(16) 
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From  the  above,  there  are  many  lists  (L,A,A, A) .  The  primary  index  for  these 
lists  is  based  on  the  knowledge  that  the  Kalman  filter  provides  the  filtered 
estimate  for  sensor  measurements  from  a  target.  Because  of  this  and  the  fact 
thaimore  than  one  target  is  considered,  a  second  index  or  subsorrpt  rs  also 
added  for  referencing  a  particular  list  L,  le  L.  .  .  The  first  subscript 

denotes  the  position  of  the  origin  with  2^  -1  in  L.  Thus  i-2  represents 

case  of  the  second  smallest  d  value  originating  from  the  reference  target 
The  range  of  the  subscript  '  i’  is  i=0,..,m  where  the  list  LQ.  contains 

origin  equal  to  the  reference  target  index.  The  second  subscript _  'j'  denotes 
the  permutation  of  the  remaining  m-1  origins  with  2.>1.  Thus  j  represents 

an  index  for  track  origin  arrangements  and  has  the  range  j=l, Examples 

of  the  track  origin  arrangements  LQj  and  for  m=2  and  t=2  are, 


L01 

=  (0,0} 

L02 

=  (0,2) 

L03 

=  (2,0} 

list  L..  the  following  s 
ij 

=  number  of 

clutter 

11 


12 


ij 


ij 


L  =  set  of  target  origins  not  in  L , 

ij  •  J 


j=l, • • 


(17) 


(18) 


Note  that  elements  z  associated  with  a  list  L.  .  have  components  with  some  o£ 
the  origins  in  L...  These  elements  have  the  d  value  associated  with  an  origin 
tfL  greater  thin  dt  ,  ie  not  in  the  m  smallest.  Therefore  all  elements  z 

which  have  the  d  value”of  a  measurement  with  target  origin  1.  in  the  m  smallest 
are  represented  by  the  event  M,.  By  adding  the  measurement  origin  to  the  d 

value,  the  event  is  defined  by, 


M0  =  U;  (d.,2  =2)<(dx  )  and  VzsD£nV£> 

2  11  mm 


11=0,1,  .  .  ,t 


(19) 


Therefore  from  equation  (19),  all  elements  in  DjnVjnM,  have  the  d  value  of 
the  measurement  with  origin  l  outside  the  m  smallest  but  contained  in  the 
element  z. 

From  the  preceding,  each  element  z  corresponding  to  the  origins  in  the  set  L.. 
has  n-m  origins  which  contain  none  or  any  one  of  the  possible  target  ongxn 
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combinations  selected  from  the  origins  in  .  Each  target  combination  is 

contained  in  the  lists  (A, A).  When  an  element  z_  has  an  origin  izK  and  Izz 
then  zsD^nV^nM^.  Conversely  when  the  origin  Iz A  and  Zfz  then  zeD^TTv^.  To 

represent  these  possibilities,  the  binary  number  defined  in  (14)  is  utilized. 
In  this  case  the  number  of  bits  is  t-m+0^  which  is  the  number  of  target 

origins  in  L_  .  Therefore  the  binary  number  is  B^Ct-m+0  ).  The  members  of 

A  are  now  obtained  from  the  direct  product  of  the  members  of  B  (t-m+d  ) 
1JK_  jc  ij 

and  L  .  Clutter  origins,  ie  zero  entries,  are  not  retained  in  A  .  The 
1J  ijK 

members  of  A^^  are  therefore  obtained  from  complementing  the  binary  number 

before  performing  the  direct  product  with  L .  These  sets  and  the  associated 

symbols  to  describe  them  are  defined  by 


A..  =  target  origins  from  the  direct  product  of  the 

x  J  _ 

members  of  L. ,  and  B  (t-m+d,.), 
ij  ij" 

6.  =  number  of  bits  with  b.=0  in  B  (t-m+b  ). 

XJK  i  ij 

A  =  the  6  target  origins  in  L. .  U  A.  .  ,  <=1, . . ,3 . 

]K  UK  IJ  11IC  ’  ’  - 


1J> 


1JK 


IJ 


(20) 


where  3..  =  2t  m+<^ij  . 
ij  J 

For  each  set  of  lists  (L  ,A  ,A  ),  there  are  many  arrangement  lists  A  to 

account  for  all  possible  permutations  of  clutter  subscripts.  These 

permutations  are  indexed  by  a  fourth  subscript  1 h r ,  ie  A.  .  The  range  of 

i  j  Kn. 

h,  which  is  represented  by  n ± ,  is  found  from  the  number  of  clutter  origins  in 

^ij  an<^  total  number  of  clutter  origins  from  n  and  the  lists 

(L. . , A . .  ,A . .  ). 
ij’  ijk’  ijk' 

From  the  definition  of  the  lists  (L,A,A, A)  and  event  M£,  the  event  space  Z 
consists  of  the  partitions 


fiijK  -  {z;  (L’A)=(Lij>AijKh)  f°r  h=l,..,TiijK  and  V  zeD..^ 


i=0 ,  .  .  ,  m ,  j=l,..,^i  and  k=1,..,3 


(21) 


where  D .  . 

ij* 


11 


ERL-0381-TR 


Figure  4  represents  diagrammatically  these  partitions  on  the  sample  space  Z. 
The  D,  V  and  M  events  on  Z  are  not  shown  because  of  the  complexity  of  the 
partitions  as  illustrated  in  equation  (21). 


From  equation  (21),  each  partition  contains  permutations  of 

subscripts.  Each  permutation  is  represented  as  a  sub-partition*  on 

Q , .  .  The  range  of  the  subscript  h  depends  on  the  D,  V  and  M  events  in  which 

1JK 

an  element  z  is  contained.  From  the  list  definitions  for  partition  Q, .  .  ,  the 

number  of  clutter  measurements  in  L. .  is  0 . .  and  the  total  number  of  clutter 

ij  iJ 

measurements  from  which  these  are  selected  is  n-t+6. .  .  Thus  the  number  of 

ljlC 

clutter  index  permutations,  and  hence  the  number  of  sub-partitions  is 

(n-t+6 . .  ) ! / (n-t+6 . .  -0  .  . ) ! . 
ijk'  ijk  ij' 

From  the  above,  the  mapping  from  Z  to  Y  by  the  function  iKz)  (see  figure  4) 
involves  a  re-arrangement  of  subscripts  and  is  a  many-to-one  mapping[10].  The 
inverse  from  each  element  ys8.  .  consists  of  ti  .  ,  branches  onto  Q.  .  .  Each 

J  IJK  IJK  1JIC 

branch  is  a  one-to-one  from  an  element  ye9..  to  an  element  z  of  a 

IJK 

sub-partition  The  partitions  or  events  on  the  ordered  space  Y  are 

defined  as 


0.,  =  (y=f  (z) ;  L=L.  .  and  Vz  } 

ijK  J  1J  IJK 

i=0 ,  .  .  ,  m ,  j=l,..,$i  and  K=l,..,a„ 


(22) 


The  union  of  9 . .  over  all  i,i  and  k  gives  the  ordered  space  Y.  The  union  of 
IJK 

0..  over  all  j  and  k  gives  the  partitions  (or  events)  0.  defined  by, 
i  J  k  r 


«i 

0.  =  U 

i 

j=l 

a .  . 

ij 

U  0  .  . 

IJK 

K—  1 

p. 

II 

o 

S 

=  {y; 

*.=1} 

1 

i=l, . . ,m 

=  (y; 

all  H^l} 

i=0 

(23) 

This  sub-partition  is  only  strictly 
include  the  equal  operator,  ie  d^<d^.  . 

correct  when  the  ordering  does  not 
In  this  paper,  the  ordering  defined 

in  equation  (8)  is  used  and,  since  the  probability  of  two  or  more  d  values 
being  equal  is  effectively  zero,  each  arrangement  is  considered  to  form  a 
sub -part it ion. 
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The  events  0^,  i=l,..,m  then  contain  all  elements  y  which  have  the  i  nearest 

measurement  originating  from  the  reference  track.  The  event  0^  contains  all 

those  elements  y  which  do  not  have  measurements  from  the  reference  track  in 
the  m  nearest. 


5.  TARGET  STATE  ESTIMATE  FOR  PDA  FILTER 

From  Section  3,  the  values  output  by  the  sensor  at  time  'k'  for  the  reference 
track  are  (Y(k) ,n(k) ,T(k-l) ) .  From  the  definition  of  the  sample  space  Y,  many 
elements  can  be  found  which  have  components  with  these  values.  Also  the 
elements  which  contain  these  values  are  distributed  throughout  the  events  or 
partitions  0  ^ , ©m  on  Y.  The  state  estimate  for  the  target  associated  with 

the  reference  track  is  then  given  by  reference  2. 


x(k|k)  =  E{x(k) |Y(k) ,n(k) ,T(k-l))} 
m 

=  ^  e.E{x(k)|0i(k),Y(k)>n(k)JT(k-l))}  (24) 

i=-l 


where  0.  =  Pr{0.(k) |Y(k) ,n(k) ,T(k-l))} 

That  is,  the  state  estimate  is  the  weighted  sum  of  the  estimates  of  the 
target's  state  conditioned  on  all  the  target  data.  The  weighting  term  is  the 

probability  that  the  ith  selected  measurement  is  from  the  target.  At  this 
stage  it  is  assumed  that  the  event  probabilities  0^  are  available.  Note  that 

the  density  associated  with  the  estimate  involves  the  sum  of  weighted  normal 
densities.  These  are  obtained  by  track  split  back  to  the  start  of  the 
track(ref .9) .  To  avoid  the  complexity  of  this  procedure,  the  following 
assumption  is  introduced. 

ASSUMPTION  NO. 7:  An  approximate  sufficient  statistic  for  the  past  target 
data  T^(k-l)  at  time  'k'  is  mean  x^(k|k-l)  and  covariance  P^(k|k-1). 

This  is  the  key  assumption  for  the  PDA  filter.  Using  this  assumption. 
Appendix  I  gives  the  derivation  of  the  target's  state  and  the  associated 
covariance.  The  reference  target's  state  is  given  by, 


x(k|k)  =  x(k|k-l)  +  Wy  (23) 

where 


m 


i=l 
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and  the  associated  covariance  is  given  by, 


m 

P(k|k)  =  B0P(k|k-l)  +  [1-PqJP1  Ck|k)  +  W[^  &iyiyTi-yyT]WT  (26) 

i=l 

The  covariance  in  equation  (26)  consists  of  the  sum  of  the  covariance  when  all 
the  measurements  are  clutter  and  the  covariance  when  the  correct  measurement 
updates  the  filter.  The  final  term  in  the  sum  reflects  the  uncertainty  of  the 
weighted  innovation  y  by  adding  its  associated  covariance. 


6.  EVENT  PROBABILITIES  FOR  THE  PDA  FILTER 

The  derivation  of  an  expression  for  the  event  probabilities  g^(see 

equation  (24))  requires  expressions  for  the  density  function  of  an  arbitrary 
y,  the  probabilities  associated  with  the  events  D,  V,  and  M,  and  the 
probability  associated  with  the  partition  of  Y  by  the  value  n(k) .  This 
section  gives  an  outline  of  the  derivation  of  these  probabilities. 

From  Assumption  No.l,  the  density  of  y  for  a  clutter  measurement  in  the 

region  V  is  1/v  .  From  Assumption  No. 7,  the  probability  density  of  a 
C  G 

measurement  from  target  ' i'  at  time  ' k '  using  data  from  time  k-1  is  normal. 

With  y  unconstrained,  ie  not  restricted  to  being  in  the  region  V  ,  the 
l  ° 

density  for  an  arbitrary  clutter  measurement  and  an  arbitrary  target 
measurement  is  then  given  by, 


fo(yi) 


W 


p(yijili=0) 

i/vc  if  y±  ^  vc 

0  otherwise 

p(yi>*i=iMT£) 

exp^y7.^"1?.^/  [(2ir)1J/2/|Sjl|]  £=l,..,t  (27) 


where  p  =  number  of  components  of  vector  y. 

The  probability  for  event  M,  which  is  presented  later,  requires  the 
probability  Prided}  for  an  unordered  d^  For  clutter  measurements,  this 

probability  is, 


GQ(d)  =  Pr{(d.,Z.=0)<d} 


=  v(vdnvc)/v 


c 


0  <  d  <Z 


(28) 
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/vT  “  l_/v 

where  v(.)  denotes  volume,  V,  is  the  multidimensional  ellipsoid  {y  S  y^d) 

a 

with  volume  v^ir^^d2^  |  S 1 2/T  (p/2+1)  where  p  is  the  number  of  elements  in  the 
measurement  vector  y.  X  is  the  largest  d  which  is  contained  within  . 

In  the  variable  'd'  when  V ,c  V  ,  the  density  for  the  reference  target,  1=1, 

d  c 

is  chi-squared  with  p  degrees  of  freedom.  The  probability  Pr{d^<d)  for  a 
target  measurement  is  then 


Gx(d)  =  PrUd.^l^dlTp 


f  .  p/2-1  -u/2  , 

u  e  du 

0  2y/2r(p/2) 


(29) 


where  the  term  is  a  scaling  factor  to  account  for  the  intersection  of  the 
shell  with  V^.  In  many  applications  this  does  not  occur  and  therefore 
A^=l.  For  series  approximations  to  the  integral  in  equation  (29)  see 
reference  8. 

For  other  targets  in  the  cluster,  ie  £>1,  no  similar  elegant  function  is 
available  for  Prld^Sd).  Therefore  it  is  assumed  that 

ASSUMPTION  NO. 8:  The  probability  function  Pr{ (d^  , £+=£)^d | T^} ,  for  Z> 1  is 
available  and  represented  by  G^(d). 

An  expression  for  G^(d)  involves  the  integral  of  the  multivariate  normal 

density  function  f^(y)  over  the  offset  region  V^.  An  approximate  expression 

for  this  probability  is  obtained  by  assuming  that  the  covariances  of  track  ' £' 

and  the  reference  track  are  related  by  ^=c  ^"S  \  where  c  is  a  scalar. 

Under  these  conditions,  the  region  can  be  made  circular  by  scaling  the  y 

coordinate  space  so  that  track  ’  has  unit  variance.  Thus  the  density  of  d/c 
is  non- central  chi  squared  with  p  degrees  of  freedom  and  non-centrality 
parameter  /t ^/c .  Series  approximations  for  this  distribution  can  be  found  in 

reference  8. 

The  event  probabilities  for  D,  V  and  M  can  now  be  given.  When  these  events 
are  assumed  to  be  independent  of  each  other  and  of  n  (see  Assumption  No  s  1 
and  3),  they  can  be  written  as 
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Pd£  =  pr{D£(k)|T£(k-l)} 

Pv£  =  pr{V£(k)|D£(k),T£(k-l)} 

Pm£  =  Pr{M£(k)|D£(k),V£(k),T£(k-l)>  1=0,.., t  (30) 


If  any  of  the  above  events  are  not  independent  of  n,  then  n  has  to  be  included 
on  the  LHS  of  the  conditioning  and  the  expression  for  the  event  probability 
then  includes  n. 

From  the  definition  of  the  above  events,  the  probabilities  for  origin  value 

£=0  ie  P,.,  P  _  and  P  A  are  unity.  The  probability  of  target  detection, 
d0  vO  mu 

P^£,  is  application  dependent  because  it  is  a  function  of  the  detection 
processing  illustrated  in  figure  2.  Therefore  the  assumption  is  made  that 

ASSUMPTION  NO . 9 :  The  event  probability  P^£  for  £=l,..,t  is  known. 

The  value  for  the  probability  of  a  measurement  from  the  target  being  in  the 
region  V  ,  ie  P  ,  is  obtained  from  the  integral  of  f  (y . )  over  the  region  V  . 

This  probability  can  be  obtained  by  series  approximations  and  from  the 

dimensions  of  V  . 

c 

The  probability  of  a  measurement  from  target  £  having  a  d  value  in  the  m 

nearest,  P  is  equal  to  G. (d)  with  the  value  d=d  (see  equations  (28),  (29) 
m£’  £  m 

and  Assumption  No.  8).  Thus  pmjj=G£ (dm)  • 

The  last  remaining  probability  required  to  obtain  a  solution  for  the  event 
probabilities  deals  with  the  partition  of  0  by  the  number  of  measurements 

n(k) .  From  Section  4,  the  elements  y  have  n  ranging  from  m+1  to  N1  .  The 
partition  caused  by  n(k)  is  expressed  in  terms  of  the  probability  of  obtaining 
'r'  clutter  measurements  in  the  region  given  T(k-l).  This  probability  is 

represented  by  $(r|T(k-l))  and 

ASSUMPTION  NO .10:  The  a  priori  probability  of  having  n  measurements  from 
clutter  is  equal  to  that  for  n-0  measurements  from  clutter, 
ie  f(n|T(k-l))=#(n-*|T(k-l))  where  0=1,.., t. 

An  alternative  to  this  assumption  is  to  adopt  the  approach  taken  in 
reference  3  where  a  Poisson  distribution  with  known  clutter  measurement 
density  is  used.  This  approach  is  not  taken  here  because  n  is  assumed  to  be 
large . 

From  the  events  defined  in  Section  4,  the  preceding  density  and  probability 
functions,  the  event  probabilities  can  be  derived  as  shown  in  Appendix  II,  and 
are  given  by 


m 


j=0 


i=0 , . . ,m 


(31) 
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The  terms  are  obtained  by  evaluating, 


r. 


b. 

1 


3  .  . 

ij 


j=l  £  eL. . 
i-  ij 


K=1 


where 


1  Pde(Pv2'Pmt)1  11  “t» 


£sA .  . 


£sA .  . 

1JK 


aijK 


(n-t+6 . .  ) ! 
1JK 

(n-t+6  .  .  -0  .  . ) ! 
1JK 


(1-P..P  „)/p,„p  „ 
d£  v£  d£  v£ 


(32) 


Equation  (32)  contains  a  density  which  is  proportional  to  the  product  of  the 
joint  density  for  the  targets  and  clutter  measurements  in  the  list  times 

the  probability  of  having  the  other  n-m  clutter  and  target  measurements 
outside  the  m  nearest.  This  joint  density  is  summed  over  all  j  origin  lists. 
Note  that  from  equation  (32),  the  time  dependent  functions  are  n,  f^  (yi), 

P,  ,  P  „  and  P  „ .  An  example  of  the  derivation  of  the  b.  terms  for  the  case 
d£  v£  m£  i 

m=2,  t=2  is  given  in  Appendix  III.  The  validity  of  these  results  can  be 

checked  by  setting  P^^O  f°r  the  terms  in  Appendix  III.  When  the  b^  terras 

are  all  multiplied  by  PV^P^2  and  common  terms  cancelled,  the  following  is 

obtained. 

b0  =  n“lPvl/Vo  +  <"‘2)(PvrPnl)/Ivc(1'PmO)1 

bl  =  fl^Yl^  ("33') 

b2  =  fl(y2) 

This  is  the  same  result  as  that  given  in  refernce  7  when  track  initiation  is 
not  included  and  when  P  =1.  When  Pd2=1  and  f2(y)=0,  the  b0  term  chan§es  to 

b0  =  (n-1)“lPvl''vc  +  (n-3)(pvrpml)/[vc<1-p„0:>' 


(34) 
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This  is  the  correct  result  since  in  this  case  the  number  of  clutter 
measurements  is  reduced  by  one. 

From  the  above  it  can  be  seen  that  if  there  are  not  ' t *  tracks  in  the  region 
V  ,  dummy  tracks  with  a  low  P^  (ie  0.001)  will  not  result  in  any  interference 

with  the  valid  tracks. 

As  can  be  seen  from  the  development  given  in  Appendix  III,  the  algorithm 
complexity  rapidly  increases  with  'm'  and  ' t ' .  No  general  expression  has  been 
found  for  deriving  the  number  of  events  from  'm'  and  't'.  Therefore  the 
events  are  written  down  from  the  origin  lists  (see  definitions  (16),  (18)  and 

(20)).  For  most  practical  applications,  t=3  and  m=3  are  reasonable  upper 
limits.  For  this  case,  when  i=0,  £=13  and  when  i>0,  £=7 .  The  intermediate 
case  of  t=2  and  m=3  has  £  values  of  4  and  3  respectively.  While  this 
complexity  exists,  the  algorithm  is  the  same  for  all  tracks  and  there  is  no 
on-line  requirement  for  the  computation  of  the  feasible  events  as  in 
reference  3,4.  The  algorithm  adapts  to  changing  clutter  conditions  via  the 
terms  n  and  P  . . 


7 .  CONCLUSIONS 

The  use  of  gates  to  both  select  sensor  measurements  and  locate  tracks  which 
are  likely  to  interfere  results  in  a  PDA  algorithm  whose  number  and  complexity 
of  computations  change  with  conditions.  This  paper  has  developed  the  theory 
for  a  PDA  algorithm  based  on  alternative  rules  for  selecting  measurements  and 
interfering  tracks.  These  rules  replace  the  requirement  for  a  gate  by  a 
nearest  neighbour  rule  where  more  than  one  neighbour  is  selected.  While  the 
algorithm  is  more  complex  than  a  PDA  algorithm  with  tracks  in  isolation,  it 
remains  unaltered  with  clutter  and  target  conditions. 

Because  tracks  are  all  updated  with  the  same  computations,  an  array  processor 
can  be  used  to  achieve  the  required  processing  power  for  the  application  of 
this  algorithm  to  real-time  tracking  on  H.F.  and  microwave  radar  systems.  The 
algorithm  presented  in  this  paper  is  to  be  extended  in  a  subsequent  paper  to 
include  non-uniform  clutter  density,  track  initiation  and  termination. 
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APPENDIX  I 

DERIVATION  OF  THE  PDA  FILTER 

An  expression  for  the  state  estimate  of  each  reference  target  requires  the 
mean  E{x|0^,Y,n,T} .  This  mean  is  considered  separately  for  the  cases  of  i=0 

and  i=l , . . ,m .  From  Assumption  No. 7,  the  Normal  density  function  applies  and 
it  is  abbreviated  in  this  appendix  as  N{mean, variance} . 

For  i=0 ,  the  mean  given  the  event  that  all  measurements  Y  are  not  from  the 
reference  track  and  past  track  data,  is 

E{x|0o,Y,n,T}  =  E{x | Ti } 

=  x(k|k-l)  (1-1) 


Note  that  in  this  case,  no  use  can  be  made  of  data  Y  and  the  current  target 
state  is  independent  of  other  targets  (see  Assumption  No. 3).  The  density 
associated  with  the  estimate  (see  Assumption  No. 7)  is 


N{x(k|k-1) ,P(k|k-l) } 


(1.2) 


th 

For  i=l , ■ . ,m,  the  mean  given  that  the  i  measurement  is  from  the  reference 
track  and  all  other  measurements  are  not,  is  as  before  based  only  on  the 
past  data  for  the  reference  track.  Thus  from  Assumption  No's  1,  2  and  4, 
and  the  Kalman  equations  (3)  to  (8),  the  mean  is  given  by; 


E{x|0.,Y,n,T} 


E(x|yi,T1) 


=  x(k|k-l)  +  Wy^ 


=  x.(k|k) 


(1.3) 


For  equation  (1.3)  the  target  index  in  equation  (3)  has  been  omitted  because 
JL=1.  The  density  associated  with  this  estimate  is  that  for  the  Kalman 
filter  and  is  represented  by 


N{x.(k|k),P1(k|k)>, 


(1.4) 


where  the  covariance  is  obtained  from  equation  (4) .  Note  that  the 
covariance  is  independent  of  the  measurement  'i'  because  it  is  conditioned 
on  this  measurement  being  from  the  target. 

From  equation  (23),  the  event  probabilities  sum  to  unity.  Therefore  with 
equations  (1.1)  and  (1.3)  substituted  into  equation  (24),  the  state  estimate 
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of  the  reference  target  is  given  by; 


x(k|k) 


m 


PQx(k|k-l) 


i=l 


e.x.(k|k) 


=  x(k|k-l)  +  Vy 


where 


(1.5) 


Now  from  equation  (24) ,  the  covariance  associated  with  this  estimate  is  given 
by; 

P(k|k)  =  f  [x-x(k|k)][x-x(k|k)]Tp(x|Y,n,T)  dx 


[x-x(k|  k)  ]  [x-x(k  |  k)  ]Tp(x  |  0_^ ,  Y,n,T)  dx 


(1.6) 


By  substituting  [ (x-x^(k | k) )  +  (x^ (k | k) -x(k  |  k)  ) ]  for  [x-x(k|k)]  in 

equation  (1.6)  and  by  using  the  density  defined  in  expression  (1.4),  the 
integral  in  equation  (1.6)  for  i=l,..,m  becomes, 

l  [x-x±(k|k)] [x-xi(k|k)]TN{xi(k|k),P1(k|k)}  dx 

+  [x.(k|k)-x(k|k)]  [x.(k|k)-x(k|k)]T  j  N{x.(k|k)  .P^klk)}  dx 
+  [x.(k|k)-x(k|k)]  |  [x-x.(k|k))]TN{x.(k|k),P1(k|k)}  dx 
+  f  [x-xi(k|k)]N{xi(k|k),P1(k|k)}  dx  [x^ (k | k) -x(k | k) ) ]T  (1-7) 


The  first  integral  in  equation  (1.7)  evaluates  to  P1(k|k),  the  second  is 
unity,  and  the  third  and  fourth  are  zero.  Similar  results  are  obtained  for 
i=0.  Therefore  equation  (1.6)  can  be  written  as. 
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m 


m 


P(k|k)  = 


B  P(k|k-l)  + 


B^P1 (k|k) 


B±x  (k|k)x  (k|k) 


i=l 


i=l 


x(k|k)x"^(k|k) 


Cl. 9) 


From  the  definition  of  x(k|k)  and  x  (k|k)  in  equations  (1.5)  and  (1.3),  The 
covariance  of  the  state  estimate  can  be  written  as, 


m 


.  X  /v  /vT  A/  a/T  ■,  t  tT 

=  B0P(k|k-i)  +  [ i-BqJp1  (k | k)  +  w[  £  i-yy 

i=l 


P(klk) 


(I. 10) 
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APPENDIX  II 

DERIVATION  OF  THE  PDA  FILTER  EVENT  PROBABILITIES 

The  event  probabilities  0^  for  the  set  of  events  defined  by  equation  (24)  can 
be  expanded  using  Baye's  rule  to  give: 

Pi  =  Pr{0i|Y,nsT} 


p(Y,0i,n|T) 


m 


(numerator) 


i=0 


(II. 1) 


The  joint  density  in  the  numerator  of  equation  (II. 1)  is  the  probability 

th 

density  for  the  ordered  values  in  Y  having  the  i  value  from  the  reference 
track  and  there  being  n  measurements  in  given  the  past  target  data  T.  From 


equation  (23),  0^  is  the  union  of  0^^ 

numerator  expands  to: 


over  all  j  and  k,  therefore  the 


5.3.. 

i  ij 


p(Y,0.,n|T)  = 


p(Y,0..K,n|T) 


(II- 2) 


j=l  K=1 


From  the  definition  of  0 ^ ^ ^  in  equation  (22),  the  density  on  the  RHS  of 

equation  (II. 2)  is  for  the  values  in  Y  with  the  origins  defined  by  L__  with  n 

measurements  in  V  and  with  target  measurements  in  the  partition  D..  . 

c  ij  *• 

The  joint  density  of  the  unordered  y  values  in  Y  with  the  origins  defined  in 
L..  and  in  the  partition  Cl  (D  0  V.  fl  M  )  is  derived  on  Z  using  the  values  in 


ij 
Y,  ie, 


ZsL. . v  Z 
ij 


Z  eL.  . 
i-  iJ 


(II. 3) 


From  the  definition  of  D_k  in  equation  (21),  the  probability  associated  with 
the  target  origins  in  A . .  and  A.,  on  Z  is 

&  &  1JK  1JIC 
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[  n 

£eA .  . 

1JK 


Pdt(Pvl'P«£)1  n 


£eA .  . 

ij* 


(II. 4) 


The  probability  of  the  remaining  n-t+6  measurements  being  from  clutter 

and  outside  the  m  nearest  is 


§  (n- 


t+d.  . 
1JK 


■hj)(1 


n  (n-t+6 . .  ) 

■Pn,0>  1J 


(II. 5) 


From  the  definition  of  the  ordered  space  which  is  illustrated  by  the  mapping 
functions  in  figure  4,  there  are  (n-t+6  )  !  /  (n-t+6..  ^-0^  )  !  branches  from  the 

ordered  space  Y  to  the  unordered  space  Z.  Thus  the  density  on  the  RHS  of 
equation  (II. 2)  is  obtained  from  the  product  of  this  number  of  branches  with 
expressions  (II. 3),  (II. 4)  and  (II. 5)  to  give, 


p(v,e1:jK,nm  -  11  VVpdtll(1-p»o><6iJ,r*ij) 

i  eL.  . 

I 


[  n 

£eA  .  . 
1JK 


W'W1  " 


It  A .  . 
1JK 


(II. 6) 


where 


a 


ijk 


c .  . 

iJK 


’(n-t+6  .  .  )  ! 

IJK 

(n-t+6  .  .  -0  .  .  )  ! 
1JK  1J 


t 


*(n-t+6ijK.-0i.)(1-PmO 


)n_t  n  p  p 

}  rd£  v£ 

1=1 


(-1~Pd£Pv!l'>/Pd£Pv£ 


The  solution  to  the  0  's  now  follows  by  substituting  equation  (II. 6)  into 

equation  (II. 2)  and  then  into  equation  (II. 1).  From  Assumption  No. 10,  $(.)  is 

constant  and  therefore  c..  cancels  to  give  the  result  in  equations  (31) 

ij  K 


and  (32). 
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APPENDIX  III 

EXAMPLE  DERIVATION  OF  THE  EVENT  PROBABILITIES 


This  appendix  gives  the  derivation  of  the  event  probabilities  for  the  case  m=2 
and  t=2.  Equation  (17)  contains  the  origin  lists  L^  and  L  .  The  remaining 

lists  L„  .  follow  from  L,  .  .  Therefore  only  the  derivation  of  the  probabilities 
2j  lj 

$q  and  3^  are  given  below. 


For  i=0 ,  the  target  origin  lists  A,  A  and  parameters  0  and  5  are  given  in 
Table  III.l.  The  values  for  a^  which  are  derived  from  these  parameters 

are  also  tabulated. 


From  Table  III.l  and  equation  (32),  the  expression  for  bQ  is  given  by 


bQ  =  n(n-l)a1a2/vs2  + 

(n-l)(n-2)[(l-Fml/Pvl)<.2  +  <l->’n2/Pv2>'*ir 'V' + 

(„-2)(r1-3Hl-P|nl/PvlKl-P[n2/Pv2)/[vs(l-Pm0)l2  + 

(n-l)  [ f 2 (y -j^)  +  f2(y2)]al/(vcPv2)  + 

(.-2)[f2(y1)  +  £2(y2)Ki-Pnll/Pvl)/[vcPv2(i-Pn0)] 


(in.i) 


25 


ERL-0381-TR 


For  i— 1 ,  the  lists  A, A,  parameters  0,6  and  values  for  a^  are  listed  in 
Table  III. 2 


TABLE  III. 2  ORIGIN  LISTS  FOR  EVENT  PROBABILITY  ^ 


A 

A 

6 

a 

IC 

L  =  {1,0} 

L  =  {2} 

0 

2 

1 

(n-1) 

0  =  1 

2 

0 

0 

(n-2) 

L  =  {1,2} 

L  =  {0} 

0  =  0 

0 

0 

0 

1 

As  for  bQ ,  from  Table  III.l  and  equation  (32),  the  expression  for  b1  is 
given  by 


b1  =  (n-Df^y^Bj/Cv^)  + 

(n-2)f1(y1)(l-Pm2/Pv2)/[vcPvl(l-Pm0)]  + 

fl(y1)f2(y2)/(PvlPv2) 

For  i=2 ,  similar  results  are  obtained  and  the  expression  for  b2  is 
b2  =  (n-l)f1(y2)«2/(vcPvl)  + 

(n-2)  f  i  (y2 )  Cl-WPv2)/IVcPv1(1“PmO)l  + 
fl(y2)f2(yl)/(PvlPv2) 


(III. 2) 


(III. 3) 


partitions  W.  .  ,  formed  from  the 
r  ljich 

tutation  of  clutter  subscripts 
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Ordered  space  Y 
from  selection 
process 
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